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Abstract. We show how to reahze a quantum interface between optical fields and 
the polarized nuclear spins in a singly charged quantum dot, which is strongly coupled 
to a high- finesse optical cavity. An effective direct coupling between cavity and nuclear 
spins is obtained by adiabatically eliminating the (far detuned) excitonic and electronic 
states. 

The requirements needed to map qubit and continuous variable states of cavity or 
traveling-wave fields to the collective nuclear spin are investigated: For cavity fields, we 
consider adiabatic passage processes to transfer the states. It is seen that a significant 
improvement in cavity lifetimes beyond present-day technology would be required for 
a quantum interface. We then turn to a scheme which couples the nuclei to the output 
field of the cavity and can tolerate significantly shorter cavity lifetimes. We show that 
the lifetimes reported in the literature and the recently achieved nuclear polarization of 
~ 90% allow both high-fidelity read-out and write-in of quantum information between 
the nuclear spins and the output field. 

We discuss the performance of the scheme and provide a convenient description of 
the dipolar dynamics of the nuclei for highly polarized spins, demonstrating that this 
process does not affect the performance of our protocol. 



PACS numbers: 03.67.Lx, 42.50.Ex, 78.67.Hc 
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1. Introduction 

An important milestone on the path to quantum computation and quantum 
communication networks is the couphng of "stationary" qubits for storage and data 
processing (usually assumed to be realized by material systems such as atoms or 
electrons) and mobile "flying" qubits for communication (typically photons) [H [2]. 
Detection and subsequent storage of information is inapplicable in quantum information 
as an unknown quantum state cannot be determined faithfully by a measurement. 
Hence the development of "light-matter interfaces" that allow the coherent write-in and 
read-out of quantum information has been the subject of intense theoretical research 
[21 m [5]. Two paths have been identified to make light efficiently couple to a single 
atomic quantum system: the use of a high-finesse cavity coupled to a single atom or 
the use of an optically thick ensemble of atoms, in whose collective state the quantum 
information is to be stored. Both have resulted in the experimental demonstration of 
such interfaces [SI [71 El [HI [ID]. Even without strong coupling a quantum interface can 
be realized by combining the probabilistic creation of entanglement between atom and 
light with teleportation. This approach has been demonstrated with trapped ions [TT] . 

For qubits realized by electron spins in quantum dots [I2l [13] such interfaces have 
yet to be realized, though in particular for self-assembled quantum dots [H], which have 
many atom-hke properties, several proposals exist to map photonic states to an electron 
in a quantum dot [131 [IS] analogy to the atomic schemes. Strong coherent coupling 
between a single quantum system and a single mode of high-Q micro- and nano-cavities 
has been demonstrated experimentally [Ml [171 [HI [19], raising the prospect of coupling 
light to the quantum dot's electronic state by adapting protocols such as [3]. 

Despite their good isolation from many environmental degrees of freedom, the 
electron-spin coherence time in today's quantum dots is limited mainly due to strong 
hyperfine coupling to lattice nuclear spins. Moreover, the capacity of such an interface is 
one qubit only, making the interfacing difficult for many-photon states of the light field as 
used in continuous variable quantum information processing. In contrast, the ensemble 
of lattice nuclear spins could provide a high- dimensional and long-lived quantum memory 

m- 

We show in the following how to couple an optical field directly to the nuclear 
spin ensemble, thus interfacing light to an exceptionally long-lived mesoscopic system 
that enables the storage and retrieval of higher-dimensional states and is amenable to 
coherent manipulation via the electron spin [21]. The system we consider is a charged 
quantum dot strongly coupled to a high-finesse optical cavity by a detuned Raman 
process, introduced in Section [21 In Section [3] we show that by adiabatically eliminating 
the trion and the electron spin different effective couplings (that can be tuned on- and 
off-resonant) between light and nuclear spins are achieved. In Sections [H and [5l we 
demonstrate that the state of the cavity field can be directly mapped to the nuclear 
spins using the methods of Landau-Zener transitions [22l [23] and stimulated Raman 
adiabatic passage (STIRAP) |2l], respectively, where the latter yields a reduction of the 
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time required for write-in. However, the drawback of this approach is that it requires 
very long cavity hfetimes. To address this problem we discuss in Section [6] at length 
an approach that was proposed in [25j (here, we discuss it for an experimentally more 
promising setup), which is robust against cavity decay: the read-out maps the nuclear 
state to the output mode of the cavity, while the write-in proceeds by deterministic 
creation of entanglement between the nuclear spins and the cavity output-mode and 
subsequent teleportation [26]. In Section [6l we give further insight into this system and 
its dynamics: we describe the full time evolution of the system, compute the read-out 
fidelity and derive the shape of the output mode-function. Moreover, we show that 
apart from mapping light states to the nuclear spins, the interaction we describe can be 
used to generate an arbitrary Gaussian state. In Section [7] we discuss different aspects 
concerning the experimental realization and the approximations used in our scheme such 
as the internal nuclear dynamics, dominated by dipolar interactions, which we model 
numerically and corrections to the first order bosonic description. 

2. The system 

We consider a self-assembled III-V quantum dot (QD) with a single conduction-band 
electron strongly coupled to a high-Q nano-cavity [see figured^)]. At zero magnetic field, 
the two electronic ground states |±l/2) (s-type conduction band states) are degenerate 
and the only dipole allowed transitions are to the trion states |±3/2) with spin -1-3/2 and 
spin —3/2 (heavy-hole valence band state) with polarized light. An external magnetic 
field Bz in 2;-direction, perpendicular to the growth (y-) direction, Zeeman splits the two 
electronic states and the trion states and leads to eigenstates |±) = -^(|l/2)±i| — 1/2)) 
and |T±) = ^(|3/2)±i |-3/2)). The states |+) ^ |T+) and |-) ^ |T_) can be coupled 
[see figure [T)d)] by horizontally polarized light, and |-|-) |T_) and |— ) -v^ |r+) by 
vertically polarized light: 

+ y(e+^^'*(l+)(T+| + |-)(T_|) + h.c.) 

+ ujc a^a + UJT+ |T+)(T+| + ujt_ \T-){T^ \ + uJeS\ (1) 

Here, h = 1, is the electron spin operator, and a are creation and annihilation 
operators of the single mode cavity field and c^c, denote the cavity and the laser 
frequency (which are vertically/horizontally polarized respectively) and fic, the Rabi 
frequencies of the cavity and the laser field, respectively. The energies of the trion states 
|T+), |T_) are ujt+ = ujT + uJh/'2, and ujt_ = ujT — oJhl'2. where ut is the energy of the trion 
(without magnetic field), the energy of the hole Zeeman splitting and = g^fibBy 
denotes the Zeeman splitting of the electronic states. The first term of the Hamiltonian 
given by ([1]) describes the coupling to the cavity field and the second term the coupling 
to a classical laser field. We assume both cavity decay and spontaneous emission rate of 
the QD to be much smaller than flc and omit both processes in ([T]). Besides the coupling 
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to optical fields, the electron spin in a QD also has a strong hyperfine interaction with 
the lattice nuclear spins, which is for s-type electrons dominated by the Fermi contact 
term 

Hm = ^i^^S- + S+A-) + AS' A', (2) 

where S^'^ are the electron spin operators and A^'^ = V . ajlf'^ are the collective 
nuclear spin operators (in a typical GaAs quantum dot, the number of Ga and As nuclei 
lies between ~ 10^-10^). The individual coupling constants aj are proportional to the 
electron wave function at site j (and the magnetic dipole moment of the jth nucleus) 
[27] and are normalized to aj = 1 . The requirement for using nuclear spins as a 
quantum memory is to initialize them in a well-defined, highly polarized state. By this 
we mean that {A') is close to its minimum value (^^)jnm (~ ~V2 for spin-1/2 nuclei) 
and define the polarization as P = (A') / {A')^^^. Due to their small magnetic moments, 
nuclear spins are almost fully mixed even at dilution-fridge temperatures and fields of 
several Tesla. Over the past years, large progress in dynamical polarization experiments 
[281 [29t [301 [31] has been reported with nuclear polarization up to 60%, recently, nuclear 
polarization > 80% has been achieved [32] . 

A convenient and intuitive description of the highly polarized nuclei with 
homogeneous coupling to the electron is provided by the Holstein Primakoff 
transformation [33], by which collective nuclear spin operators A^'^ can be mapped 

1 . h m u ^^A Az , 1 fhih 



to the bosonic operators b,b\ associating A" VWy'^ ~ ^ Jj (^^^ ~ t)- 

Assuming high polarization, the electron spin couples to a bosonic "spin wave" described 
by A~ = -^b and A' = -^{b^b — N/2) by a Jaynes-Cummings-like interaction 

= |(&^5- + Sn) + -^S' (b^b - 1) , (3) 



with gn = A/yN. The initial state of the nuclear spins is represented by a collection 
of bosonic modes, with b in the vacuum state. One can generalize this description 




Figure 1. (a) A singly charged QD coupled to a high-Q optical cavity, (b) Level 
scheme of the QD. Optical and hyperfine driven transitions. 
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to the case of inhomogeneous coupling to the electron = A^y^■ af) and obtains 
an identical description in 0th order in {^^^ = (l~-P)/2 [34j. Corrections to this 
description arising from inhomogeneous coupling and not fully polarized nuclear spins 
will be discussed briefly in Section fL2\ and a detailed discussion can be found in [35]. It 
should be noted that the scheme we present does not require the bosonic description and 
could also be discussed directly in terms of the collective spin operators. The Fock basis 
would be replaced by (A"'")"' || ... J.) and errors due to the inhomogeneity would have to 
be treated along the lines of [20] and [25]. The bosonic picture, however, allows a much 
more transparent treatment of the corrections to the ideal case, emphasizes the relation 
to quantum optical schemes, and gives access to the Gaussian toolbox of entanglement 
criteria and transformations. 



3. Effective coupling between nuclei and cavity 

Our aim is to obtain from H = Hopt + Hm a direct coupling between nuclear spins and 
light. The Hamiltonian H describes a complicated coupled dynamics of cavity, nuclei 
and quantum dot. Instead of making use of the full Hamiltonian (and deriving the 
desired mapping, e.g., in the framework of optimal control theory) we aim for a simpler, 
more transparent approach. Closely following [25], we adiabatically ehminate [36] the 
trion and the electronic spin degree of freedom, which leads to a Hamiltonian Hei that 
describes a direct coupling between nuclear spins and light. As explained later, this 
can be achieved if the couplings (the Rabi frequency of the laser/cavity, the hyperfine 
coupling, respectively) are much weaker than the detunings to the corresponding 
transition: 

A^^ > Q,,r]cv^, (4a) 
^A^^ cDe > fi/,^^cV^, (46) 

UJe > gnV^- (4c) 

Here, A^^ = ujt — uji ± uJh/2 + uJe/2 = A' ± uJh/2 + a}e/2 with A' = ujt — uji, n is the 
number of cavity photons and m the number of nuclear excitations. Note that typically 
Ue < A'rp^ such that condition (Haj) becomes redundant. In addition to (l¥aj) -(l4cl). we 
choose the adjustable parameters such that all first order and second order processes 
described by H are off- resonant, but the (third order) process in which a photon is 
scattered from the laser into the cavity while a nuclear spin is flipped down (and its 
converse) is resonant. This leads to the desired effective interaction. 

The idea of adiabatic elimination is to perturbatively approximate a given 
Hamiltonian by removing a subspace from the description that is populated only with 
a very low probability due to chosen initial conditions and detunings or fast decay 
processes. If initially unpopulated states (in our case the trion state \X) and the 
electronic spin-up state |t)) are only weakly coupled to the initially occupied states, 
they remain essentially unpopulated during the time evolution of the system and can be 
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eliminated from the description. The higher order transitions via the eliminated levels 
appear as additional energy shifts and couplings in the effective Hamiltonian on the 
lower-dimensional subspace. 

The starting point is the Hamiltonian H = i/opt + -f^hf given by and ([2]). 
In order to get a time-independent Hamiltonian, we go to a frame rotating with 
W = exp [—iujit{a^a + |T+)(T+| -|- |T_)(T„|)] and obtain: 

H' = Y l-X^+l + «^ l+X^-l + h-c.) + y (|+)(T+| + |-Xr_|) + h.c.) 
+ 6a^a + At^ \T+){T+\ + At_ |r_)(r_| + uo,S' + H^f, (5) 

where At± = ujt^ — ui and 6 = Uc — uji. 

Choosing the cavity and laser frequencies, Uc and ui, far detuned from the exciton 
transition and the splitting of the electronic states much larger than the hyperfine 
coupling gn, such that conditions fHflj) - fj¥cj) are fulfilled, we can adiabatically eliminate 
the states \T±) and |+). A detailed derivation of the adiabatic elimination can be found 
in [25]. It yields a Hamiltonian, that describes an effective coupling between light and 
nuclear spins 

Ha = caV- M + h.c.) + ^ {aA- + h.c.) + LJ^a^a - -6A' 

^A+A'+Tnl, (6) 



4(1), 



e 



where the energy of the photons ijJ\ = 6 — — V -^^^^ and the energy of the nuclear 
spin excitations ^ ~ 4^^- By we denote the nonlinear terms 

T,i = ^A+MM- + ^ba^aA^A- + ^-a^a^aa, (7) 

which are small (||Tni|| -C f^^) in the situation we consider {6 ^ Vlc,gn/^z ~ 
VLi/A'rp^rp_ <^ 1) and neglected in the following. We also neglect the nuclear Zeeman 
term which is of order 10~^ smaller than the Zeeman energy of the electron. In the 
bosonic description of the nuclear spins that we introduced in ([3]) the Hamiltonian given 
by (P) reads 

i^eff = giio-b^ + h.c.) + g2ia'b + h.c.) + cjia^a + iU2b^b, (8) 
with coupling strengths gi and g2 given by 

_ ^c^ign _ ^c^Wn ff.^ 

The energy of the nuclear spin excitations can now be written as ti;2 = — 2^ — The 
first term in the Hamiltonian is a beamsplitter type interaction ~ (a6^-|-h.c.) whereas the 
second term is a two- mode squeezing type interaction ~ {ab + h.c). Both interactions 
can be made dominant by choosing the resonance condition to be either ui = 002 or 
uji = —UJ2. This will be discussed in detail in the following and illustrated numerically. 
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First, we validate the adiabatic elimination by a numerical simulation which 
compares the evolution of states \E'2o (where the first subscript indicates the number of 
photons and the second the number of nuclear excitations) [under the condition ui = U2, 
see figure [2^)] and \l/oo [under the condition ui = —uj2, see figure [2)d)] under the full 
Hamiltonian given by ([5]) to the evolution under the eliminated Hamiltonian given by 
(IH]). The solid lines show the evolution under the full Hamiltonian H', the dashed lines 
under the eliminated Hamiltonian iJgfj and we find that H' is well approximated by 
Hes, and that the nonlinear terms can indeed be neglected. 

For the simulation, we choose the parameters as follows: we assume a hole g-factor 
gh = —0.31 and an electron g-factor ge = 0.48 [37J; the number of nuclei = 10^, 
the hyperfine coupling constant A = lOOfieV, the laser and cavity Rabi frequency 
Qc = = the detuning of the trion A' = lOOOfieV, the effective Zeeman 

splitting of the electronic states Ue = 13.9fieV (the magnetic field in x-direction is 
4T) and the Zeeman splitting of the hole = —71.8^eV. With these parameters, 
the conditions given by (Haj) -(l46j) are fulfilled and values gi = 2.1 ■ 10~^/ieV and 
g2 = 1.9 • 10^'^/ieV are obtained. We assume full nuclear (spin-down) polarization 
and use the bosonic description. 

As already mentioned, two distinct resonance conditions are chosen in figures [2] a) 
and b), leading to different dynamics of the system: 

For resonant exchange of excitations between the two systems, we choose ivi = uj2, 
where the tuning can be done by changing 6 = ujc — uJi- Then i^cfr describes a 
beamsplitter-like coupling of the modes a and h and the effective interaction is described 
by 

^^bs = gi{ab^ + h.c.) + ujio^a + uj2b%. (10) 

Processes in which absorption (or emission) of a cavity photon is accompanied by 
a nuclear spin excitation are resonant, whereas the squeezing interaction given by 
g2{cL^h'^ + cib) is off-resonant. This can be seen going to a frame rotating with uji. g2 is 
rotating with 2uji and as 2uji ^ g2, the squeezing type interaction is off resonant. 

Tuning the energies such that Ui = —UJ2, the creation of a nuclear spin excitation is 
accompanied by scattering of a laser photon into the cavity, i.e. the effective coupling 
becomes g2{a''b'' +ab) and the beamsplitter-type interaction gilab'' + a'^b) is off-resonant. 
The driving laser now facilitates the joint creation (or annihilation) of a spin excitation 
and a cavity photon, realizing a two-mode squeezing Hamiltonian 

^^sq = 92{a''b^ + ab) + uJiO^a + uj2b^b. (11) 

The plots in figures [2^) and b) illustrate that the dynamics of the system can indeed 
be approximated by (|TOl) and ( fTTl) . To simulate the beamsplitter type coupling given 
by (fTOj) . we choose uJi = uj2 and let the two-photon Fock state '?/'2o evolve under 
the Hamiltonian given by ([8]) [see dashed lines in figure [2] a)]. Almost perfect Rabi- 
oscillations can be seen between the two-photon Fock state %Ij2q and the state with two 
nuclear spin excitations showing that g2{<i^b'^ -|-h.c.) in ([8]) can indeed be neglected. 
To simulate the squeezing-type interaction we choose oji = —002 and study the evolution 
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(a) (b) 

Figure 2. a) Evolution of the two- photon Fock state ip20 under the full Hamiltonian 
H' (solid lines) and the eliminated Hamiltonian Hcs (dashed lines) tuning the energies 
such that LUi — lu2 (beamsplitter- type interaction) . b) Evolution of the state V'oo under 
the full Hamiltonian H' (solid lines) and the eliminated Hamiltonian iJoff (dashed blue 
lines) tuning the energies such that coi = —uj2 (squeezing interaction). The dotted black 
lines show the evolution under exact two-mode squeezing (up to n = 3). 



of the state ipoo under the Hamiltonian given by (IHl) [see dashed blue lines in figure [2] b)]. 
It can be seen that the state ipoo evolves into the states ipn, ip22 and ip^^, with coupling 
strengths g2^/n, depending on the number of excitations n. We have thus shown, that 
in this case, the beamsplitter-type interaction can indeed be neglected. For simplicity, 
we restricted the number of photons and nuclear excitations to 3 in our simulation, such 
that states '?/'44 and higher excitation states do not occur and the evolution of the states 
il)22 and ^/'sa does only correspond to its evolution in a space with higher excitation 
numbers at very short times. This can be seen comparing the evolution to the exact 

two-mode squeezing which generates the state ^1 — tanh^ {.g2t) tanh" {g2t) \nn) 
for which the populations up to n = 3 are plotted in figure [2] b) (dotted black lines). 

4. Landau-Zener transitions 

To map the state of the cavity to the nuclear spins, we take advantage of a formal analogy 
between the linear two-mode interaction given by (fTOi) in the Heisenberg picture and 
the Landau-Zener problem [22l[23] . In the conventional Landau-Zener problem, initially 
uncoupled Hamiltonian eigenstates of a two-level system interact at an avoided crossing. 
This interaction is achieved by slowly changing an external parameter such that the 
level separation is a linear function of time. If the system starts in the ground state, 
the probability of finding it in the excited state is given by the Landau-Zener formula 
[221 [23]. 

Here, we invoke this idea in the Heisenberg picture to achieve a mapping of the 
photon annihilation operator a to the collective nuclear spin operator in the bosonic 
approximation fe, i.e., a ^ h (the bosonic operators a and h are initially uncoupled). In 
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the following we show that our system can be transformed to a system which corresponds 
to the standard Landau-Zener problem. 

In the Heisenberg picture, the linear two-mode interaction between the cavity mode 
and the nuclear spins in the quantum dot, given by (fTOl) . is described by a set of coupled 
differential equations for the mode operators: 

The available control parameters used to effect the change of ui — 102 are the laser Rabi 
frequency f2/ and the laser frequency uJi. We consider a linear time dependence of 

— ijJ2 = Pt (13) 

for simplicity while the coupling gi is constant. 

Denoting by «((")) the operator ua-00 +vb^oo (a normalized linear combination of 
the purely photonic a_oo and purely nuclear 6_oo) and by «*((")) = its image 

under time evolution, the two operator equations of (fT2l) can be combined into 

d , (u 




which can be written completely in terms of the mode function: 

dt\vtj \ gi 002 J \vtj ^ ' 

This is the same kind of coupled differential equation that is encountered in the Landau- 
Zener problem [221 [23] ■ Following the calculation as done in [221 [23] we find that 

\v(t = oo)|2 = K^Ui{oo)Ul{oo) = 1 - e^, (15) 
where e = e~'"'^'' and 



lim f/i(t) = -K ^,;^"" ^. e-^^-e'^*'(v^t)^^-, (16) 

2 

with 72 = ^ and the constant K = ^T^exp (— ^) (for a detailed discussion see 



Appendix B). We thus find that the cavity mode operator a is mapped to 

a+00 = Vl - e^b^oo + ea_oo, (17) 

a dominantly nuclear operator for (3 small enough so that e is small, which effectively 
means for large enough times, as uji — uj2 = Pt. 



4-1- Quality of the mapping for Fock and coherent states 

In the following we will consider the quality of the mapping within the model given by 
([T|) and ([31), other imperfections will be discussed in Section [3 To evaluate the quality 
of the mapping, we return to the Schrodinger picture. The mapping of a n photon Fock 
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Figure 3. The (solid) red line indicates the fidelity of the mapping of an amplified 
coherent state of light to the nuclear spins vs the losses of the mapping e. The fidelity 
is, even for large losses, higher than the fidelity F^i that can be achieved by classical 
means, indicated by the (dotted) black line. 



state of the cavity to the nuclei leads to a mixture of Fock states with photon numbers 
< n. The fidelity with which a n-photon Fock state is mapped, is given by 



FF{n) = {n\ tic ^ 




n 

|0X0| 













\n) 



:i-e 



2\n 



where the cavity mode is traced out. In a next step, we want to know which fidelities 
can be achieved for superpositions of number states, e.g., for coherent states. 

Coherent states are representatives of the family of Gaussian states, which play an 
important role in quantum optics and quantum information processing. A Gaussian 
state is fully characterized by its first and second moments (7,(i), where 7 is the 



state's covariance matrix and d its displacement (see Appendix A ). Since the dynamics 
generated by (fTS!) is Gaussian, the mapping can be fully characterized in terms of 
covariance matrices. The mapping of a Gaussian state of light onto the nuclear spins 
corresponds to 



(7c, 4 



map 



( (1 - e^)7c + e^7„s , Vl - e^dc + edns) 

V / 1 



(19) 



where (7c, dc)c and (7^^, dns)ns describe the states of cavity and nuclear spins, 
respectively. For a coherent state mapped to the nuclei, this corresponds to the map 

'l,VT-e^a] . (20) 



The fidelity of the mapping is given by [38 

F, = \{a I Vl^a)\^ 



exp 



(l-v^ 





2" 


j a 









(21) 



The minimal goal of a quantum interface is to achieve a better fidelity than can be 
achieved by classical means. As proved in [391 SO] the classical benchmark fidelity of 
coherent states distributed in phase space according to p{a) = ^ exp (— A|ap) is given 
by Fmax = Averaging Fc over all possible coherent input states with a Gaussian 
distribution the fidelity reads 

d^ap(a)F^ = ^ , . (22) 
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For a flat distribution with A ^ large photon numbers that lead to high losses are 
dominant and therefore Fc —>■ 0. 

A way to improve the average fidelity is to amplify the coherent state either at the 
write-in or the read-out stage, thus compensating losses due to the imperfect mapping. 
Optimal phase insensitive amplification would map (7, d) — > (^^7 + (/t^ — 1)1, i^d) [41j . 
Choosing n such that k = -^7^=5, amplification and subsequent mapping can be written 
as 

(1, a), ^ 1, ^ (.(1 - e^) + e-) 1, a)^^ . (23) 

where subscript c refers to the cavity and ns to the nuclear spins. The fidelity of 
the mapping of the amplified state p^, calculated using the relations for transition 
amplitudes of Gaussian states in [12], is given by 

Fco = («| Pm(a) |«) = det {^-^-Y^^ ^'^ ^-id^s-d.Vh.^^^s)(dr.s-d.) 

/K(2-e2)+e2 \-'/' 1 , , 

A plot of Fco is shown in figure O The fidelity of the mapping of the amplified 
coherent state is always higher than the classical benchmark fidelity = \ for a fiat 
distribution i.e., the quantum interface shows high performance even for large losses e. 



SPDC M 




gEoF 




(a) 

Figure 4. a) One part (Mi) of a two-mode squeezed state, arising from a Spontaneous 
Parametric Down Conversion source is coupled into the cavity and mapped to the 
nuclear spins of the quantum dot. Thereby, the other part (M-i) gets entangled with 
the nuclear spins, b) Plot of the Gaussian entanglement of formation (gEoF) for 
squeezing parameter r = 1 vs mapping error e. The entanglement of the nuclear spins 
with part of a two-mode squeezed state (M'i) is a decreasing function for increasing 
mapping error e. 
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4-2. Storage of an entangled state 

Up to now we have shown that it is possible to transfer Fock and coherent states of 
hght onto the nuclear spin memory. However, the ultimate test for a quantum memory 
is whether it is capable of faithfully storing part of an entangled quantum system. As 
an example for an entangled light state we consider a two-mode squeezed state where 
one of its light modes, Mi, is coupled into the cavity and mapped onto the nuclear spins 
of the quantum dot [see figure |4|l)]. To see how well the entanglement is preserved, 
we compute the entanglement between the nuclear spins and the light mode M2 using 
Gaussian entanglement of formation [13]. We find it to be a monotonically decreasing 
function of the mapping error e G [0, 1]. As seen from figure Hb) the nuclear spins of the 
quantum dot are entangled with the light mode M2. This allows a remote access to the 
memory via teleportation, required for e.g., quantum repeaters. 

4-3. Mapping time 

A timescale for the mapping can directly be found considering Hamiltonian ifbs where 
the parameters are chosen such that uJi = UJ2 and set to zero in a rotating frame: acting 
for a time t = -r j g (which is for the parameters used in Section [2] on the order of 6/is) 
it maps a ^ b and b ^ a thus realizing a swap gate between cavity and nuclear spins. 
This setting however, would in contrast to the adiabatic methods discussed in this and 
the following Sections, be sensitive to timing errors because letting ifbs act for too long 
would reverse the mapping. 

5. STIRAP 

Mapping the state of the cavity to the nuclear spins is also possible considering a system 
where only the trion states are adiabatically eliminated and elimination of the electronic 
states is not required to achieve the desired interaction. We show that, with this system, 
the process of storing a state of light to the nuclear spins can be achieved by the well- 
known technique of Stimulated Raman Adiabatic Passage (STIRAP) [2l] , which has been 
studied for multilevel systems [15] and has been demonstrated in several experiments 
|24j . This scheme allows us to coherently transfer population between two suitable 
quantum states via a so-called counterintuitive sequence of coherent light pulses that 
drive transitions of a lambda or a multilevel system. It has some advantages over the 
Landau-Zener method as the choice of control parameters is easier and less constraints 
have to be fulfilled as we do not eliminate the electronic states, which allows for faster 
mapping times compared to the Landau-Zener method. 

Note, that the main source of error, the decay of the cavity, is not considered 
here. Up to now, the experimentally achieved cavity decay rate 7 of a photonic crystal 
microcavity that couples to a quantum dot is of the order of 7 ~ 10^° ^ [45j. However, 
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(a) (b) 



Figure 5. a) Schematic view of tlie time dependence of the coupling ga and the 
constant hyperfine couphng g„. b) Schematic view of the level scheme of the system 
with initially two photons. The hyperfine coupling gn is "always on", whereas ga{t) is 
an increasing function for which ga{T) 3> gn- 



we do propose this scheme here, as cavity decay rates might improve and the scheme 
might also be used in a different setup. 

For the system proposed in Section [2] the STIRAP method is not as straight forward 
as for the setup we investigated in [25] in Section V. The reason for this is, that after 
ehmination of the trion states in the system used so far, there are two different couplings: 
9a~{.S^o) + h.c.) and ga+^S^a"^ + h.c), with ga± = Qc^i/AAt^, where the first one has 
to be made off- resonant: ga- ^ ga+, which means that Uh sets an upper hmit to the 
coupling ga+ (as the condition for the adiabatic elimination is A^^^ ^ Qc,^i)- Therefore, 
we study the STIRAP scheme for the system investigated in [25], where only the coupling 
oc (S~a^ + h.c.) is present. 

In [25], we study a singly charged QD where the electronic states are Zeeman 
split by an external magnetic field in growth/z-direction (Faraday geometry). The 
electronic state |t) is coupled to the trion state |ft) (with angular momentum +3/2) by 
circularly polarized light and the electronic state ||) is coupled to the trion state 
|Jj.) (with angular momentum —3/2) with cr^-polarized light. These transitions can 
be stimulated by a cr "'"-polarized cavity field and a o"~-polarized classical laser field, 
respectively. The trion states are mixed with a resonant microwave field, whereas 
the electronic eigenstates are unchanged as they are far detuned from the microwave 
frequency and are now both coupled to the new trion eig enstates |T^) = l/\/2(|t) - |^)) 
and |Tjj,) = l/\/2{\^) + |4)), and form a double A system (see [25] for a figure). 

In a frame rotating with the laser frequency the Hamiltonian reads 

^ = ^ (a^ \m\ - limi + h.c.) + ^(ITXT,| + |T)(T,| + h.c.) 

+ + A+ |T^)(r^| + A_ |r^XTj^| + Co^S' + H^, (25) 

where 6' = Uc — uji — oJmw and A-t = tujj, — cu^ ± Vtmw Now, we derive the Hamiltonian 
where only the trion has been eliminated. If 

A±>fi^,^]cv^ (26) 
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holds, the trion can be adiabatically ehminated. This leads to the Hamiltonian 



H, 



el 



9a {S^a 



+ h.c.) + gn {S+b + h.c.) + ^S'b^b + S'a^a + CueS' 



+ 



2A_ 



+ 



2A4 



2A_ 2A 



+ 



I TXT I 



where the coupling 



9a 



A, A_ 



(27) 



(2^ 



We thus arrive, in addition to iJhf ? at an effective Jaynes-Cummings-like coupling of the 
two electronic spin states to the cavity mode governed by 

gaiS+a + h.c.), (29) 

i.e., the absorption of a cavity photon goes along with an upward flip of the electron 
spin (and the emission of a photon into the laser mode) and vice versa. 

Next we present the STIRAP scheme: the couplings of the electronic states ||) 
and II) are given by the optical fields (ga) and the hyperfine coupling (^f^). The 
Hamiltonian describing the system is given by (1271) . It is blockdiagonal H' = 
where m denotes the initial photon number. The (2m + l)-dimensional Hamiltonian 
Hm, describing the evolution of the "m-excitation subspace" can be written in the Fock 
basis {|m, |, 0) , |m — 1, |, 0) , |m — 1, |, 1) , ...}, where the first number m — k represents 
the Fock state of the cavity, [ / ] denotes the electron spin down/up state and k the 
excitation number of the nuclear spins. In this basis, reads: 



Ago 


9aVm 




Aeo 





9nVl 

















9n\/^ 

Agi 
ga^/m - 









9a^m - 

9nV^ 







9nV2 

Ag2 



V 



(30) 



with Acf. 



6' + 



[m-k)- -^k 



uJe/2, and A^, 



+ 5'{m 



k y ' 2A_ 2A+ J 2N"' ^ei ^, "--^ 2A_ 

l-l) + ^l + uoj2ioike {0,l,2..m} and / G {0, 1, 2..m - 1}. In the following, we will 
denote the states with electron spin down as "ground states", \Gk) = \m — k, J., k) and 
the "excited states" by \Ei) = |m — / — 1, t, /). The optical fields couple states Gk and 
El with k = I whereas states with k = I + 1 are coupled by the hyperfine coupling (see 
figure lib). 

We will show in the following, that by slowly increasing the laser Rabi frequency 
and thus changing ga{t) such that ga(T) ^ gn, at the final time T (see figure 15^), an 
initial state I'?/', |,0) with no nuclear spin excitations in the quantum dot and a state 

in the cavity, evolves under the adiabatic change of H' to a state where the cavity 
is empty and its state has been mapped to the nuclear spins: 



t=o 



(31) 
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t [lo-S] 

Figure 6. Plot of the fidelity of a variety of states vs time during the adiabatic 
evolution. |1), |2), |3) denote the one, five, ten photon Fock states, respectively. \a) 
denotes the coherent state with average photon number 5. The total evolution time is 
chosen to be T = 5 • 10^'^s. 



for T — > oo. A prerequisite for the mapping is, that the ground states of the Hamiltonian 
are all degenerate within each "m-excitation"-subspace so that we can keep track of the 
phases of the individual eigenstates. This can be done by choosing the parameters such 
that Agj. does not depend on /c, which is fulfilled for (^H- 2a^+ 2a7 ~ ~2N~^^i'^ ^'^ ^'dX 
Acj, = —i^rn — uOel'^ 0- Hence, the phases 0^ = ^Gm t of the individual eigenstates 
which the system acquires during the time evolution (for perfect adiabaticity) are known 
and can be corrected, e.g., by applying a magnetic field — for a time t = ^/If^^l^ 
after the state transfer to the nuclei. Here gx and fiK denote the nuclear g-factor and 
the nuclear magnetic moment, respectively. 

5.1. Numerical integration of the Schrddinger equation 

To study the quality of the mapping of a state of the cavity to the nuclear spins, we 
numerically integrate the Schrodinger Equation given by 

4 i^w) = H\t) \m) , (32) 

where H' is the Hamiltonian given by fl27|) . The simulation computes \il){t + Dt)) = 
g-i//'(t)Dt i^j^^-j^ in T/Dt steps from t = to t = T. We assume, that the change 
of ga{t) is quadratic in time, ensuring an initially slow and a finally fast increase of 
ga{t) and ga{T) ^ gn{T), ga = ^Og-n^- Parameters are chosen as follows: we assume 
a hole g-factor gh = 2.2 and an electron g-factor ge = 0.48; the number of nuclei 
= 10^, the hyperfine coupling constant A = 10/ieV, the laser and cavity Rabi 
frequency = = 13/ieV, the detuning of the trion cujj, — cu; = 103/ieV, the effective 
Zeeman splitting (Dg = O/xeV and the microwave Rabi frequency flmw = 50fieV. The 

I It can be proven by induction, that i?„i has an eigenvalue Eq^ = Ag^Vm, i.e. that 
det {Hm — Eg^I) ~ forall m G {0, 1...00} (where m is the initial photon number). 
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fidelity of the mapping we are interested in is given by the overlap of the numerically 
evolved state p{t) = |'?/'si)(V'si| and the ideal output 

F= (V^idlp(t)lV^id). (33) 

To achieve a fidelity close to one, the total evolution time is chosen to be T = 5 • 10~'^s. 
Figure [6] shows the fidelity plotted versus time for different kind of states that will be 
discussed in the following, illustrating the different aspects of mapping. 

The one photon Fock state |1) is mapped in T = 5 ■ 10~^s with a fidelity of 
F ^ 0.99 to the nuclear spins. To see that not only population but also relative 
phases are properly mapped we have simulated an approximately coherent state |a) = 
exp ^— Yl'kLo^ 1^) with average photon number |ap = 5 and find a mapping 
fidelity of F f» 0.96. Here the known phases 0^ have been compensated. 



5.2. Error processes 

The main error processes that lead to imperfections of the fidelity are the "always-on" 
character of the hyperfine coupling, the nonadiabaticity due to finite times, non-perfect 
polarization of the nuclei and the decay of the cavity. These processes will be studied 
in the following. 

The fact that the hyperfine coupling is "always-on" leads to an "error" that is 
intrinsic to our system. Different to conventional STIRAP that uses overlapping light 
pulses, we propose to adiabatically increase the coupling ga{t) so that ga(T) ^ gn and 
therefore the mapping is imperfect as gn is constant and can not be "switched off'. 
Treating the coupling small perturbation in first order perturbation theory at 

t = T, the fidelity is found to be 

where |^°) is the ideal output state for which gn = (at t=T), and |^) = |^°) + l^"*^) + .. 
is the unnormalized eigenstate of Hm- 



Another error arises from the non-adiabaticity due to finite times of realistic 
processes. For a quantitative estimate of the time T that is needed for adiabatic passage 
to occur we use the well-known adiabatic theorem [16] and numerically compute the 
minimum time T fulfilling 

[SEmM£^%^S.Mmi<s.yt.lO,n (35) 

I I 

The left hand side of fl35|l corresponds to the probability to find the system in an excited 
state \E]^) different from (the (purely nuclear) eigenstates to the eigenvalue E = 
of H' - 0„Fg™1) and the fidelity decreases with 6^. For 1/100 < 6^ < 1/10 the 
minimum time fulfilling Equation (135|) for the mapping of one photon is in the range of 
Ifis > T > 0.11/is. 
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To get an accurate description of the errors arising from non-adiabaticity we use a 
perturbative approach to treat nonadiabatic corrections and compute the phases arising 
from non-adiabaticity [17j. As the Hamiltonian H' = @^Hm is blockdiagonal, states 
with different initial photon numbers m do not couple so that we can treat every "m- 
photon" subspace separately. Moreover, we can use nondegenerate perturbation theory 
as the groundstates \Gk) that are degenerate within each subspace are not coupled: 
{Gk'\ -^H \Gk) ■ Supposing to be at t = in one of the groundstates |0Gfe(O)) of 
H'^ = Hm — Eg^I slowly varying in time, the first order correction of the energy 
eigenvalue Eq = of H is given by 

l^k ' 

The phases (j)!^ = Jq dt which the system acquires can be found by numerical 
integration of E^^. For T = 5 ■ 10~^s and initial photon number m = 1, 0^ = —1.4 ■ 10~^ 
and for m = 2, 0^ = 0.004, respectively. Thus, as expected, the errors arising from 
non-adiabaticity are small for sufficiently long times T. 

6. Quantum Interface in the bad cavity limit 

In the previous Sections, we have shown that a quantum interface can be achieved via 
direct mapping of the cavity field to the nuclear spins of the QD. But we have also seen 
that the cavity lifetimes required for high-fidelity storage are much larger than what 
is today's state of the art, i.e., as gi^2 ^ l/Tcavity, we are, compared to the effective 
coupling, in the "bad cavity limit". A second problem with this approach is that the 
quantum information we want to map to the nuclei has to be coupled into a high-Q 
cavity. This is notoriously difficult although theoretical proposals exist that should 
avoid refiection completely. Both problems can be circumvented employing ideas similar 
to [m Ho] by using the two-mode squeezing Hamiltonian Hsq [see ( ITT]) ] (note that we 
now return to the system proposed in Section [2] for the rest of the paper). As discussed 
in [25] and elaborated in more detail below, it is possible to create entanglement between 
nuclei and the traveling-wave output field of the cavity. Then, quantum teleportation 




Figure 7. Quantum teleportation can be used to write the state of a traveling-wave 
light field onto the nuclei. 
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Figure 8. Nuclear spins of quantum dots in two distant cavities can be entangled by 
interfering the travelling wave output fields of the two cavities at a beamsplitter and 
measuring. 



can be used to write the state of another travehng-wave hght field onto the nuclei (figure 
[7]) This approach gives an active role to cavity decay in the interface and can tolerate 
a bad effective cavity as long as strong coupling is achieved in ([5]). Moreover, it does 
not require to couple the quantum information into the cavity. Similarly ifbs [ (flOll 
enables read-out, by writing the state of the nuclei to the output field of the cavity. The 
entanglement between nuclear spins and output field can moreover be used to entangle 
nuclear spins in two distant cavities by interfering the output light of the cavities at a 
beamsplitter (figure [8]). 

6.1. Entangling nuclei with the output field 

The Hamiltonian of the nuclear spin-cavity system tuned to the squeezing interaction 
f fTTj) and coupled to the environment is given by 



H = g2{a^b^ + ah) +ia J ^^4u + ^-c + J ujclc^duj, 



(37) 



where c^; are the annihilation operators of the bath and 7 the cavity decay constant. We 
have specialized (fTT!) to the case ui = —002 and transformed to an interaction picture 
[JJ] with Hq = uJi{a'^a — h'^h) + uJi ^ cl^c^duj and performed the rotating-wave and Markov 
approximations in the description of the cavity decay [51]. The quantum Langevin 
equations of cavity and nuclear operators read 

d{t) = - ig2 b\t) - |a(t) - ^c,^{t) (38a) 
h{t) = -ig2a\t). (386) 

Here, Qn describes the vacuum noise coupled into the cavity and satisfies [cin(t), 4n(^')] = 
5{t — t'). The solutions of (138 ah and ( 138 fel) are given (for t > 0) by 

a{t) = p-(t)a(0) + q{tp{0) + ^ j P-{t - T)c;^{T)dT (39a) 

Jo 

h{t) = g(t)at(0) + p+(t)6(0) + ^ f q{t - r)cl{T)dr {39b) 

Jo 

§ This maps the state up to a random (but known) displacement. It can be undone using Hhs, where 
the cavity is pumped with strong coherent light for a short time |50j . 

II As was already the case in (fTT|) all optical operators are also taken in a fram rotating with the laser 
frequency . 
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where 



P± 



4*7 



cosh (i^t) ± — sinh (z/t) 



i—e 4 ginh i/t, 



with 



+ 9l 



(40) 
(41) 

(42) 



While ( IM and (JHIS) describe a non-unitary time-evolution of the open cavity- 
nuclei system, the overall dynamics of system plus surrounding free field given by the 
Hamiltonian in (|37l) is unitary. Moreover, it is Gaussian (see Appendix A ), since all 
involved Hamiltonians are quadratic. Since all initial states are Gaussian (vacuum), the 
joint state of cavity, nuclei, and output fields is a pure Gaussian state at all times as 
well. This simplifies the analysis of the dynamics and in particular the entanglement 
properties significantly: The covariance matrix [defined by (11. 4p in Appendix A of the 
system allows us to determine the entanglement of one part of the system with another 
one. In particular, we are interested in the entanglement properties of the nuclei with 
the output field. 

The covariance matrix Fns-c-o of the pure Gaussian state of nuclear spins, cavity 
and output field and thus the covariance matrix Fns-o of the reduced nuclei-output field 
system can be found by analyzing the covariance matrix of the cavity-nuclei system 
r 

ns— c- 

The elements (X) of the covariance matrix Fns-c can be calculated by solving the 
Lindblad equation evaluated for the expectation values (X) 



d_ 

dt 



{X) = i {[H,^,X]) + I {{2a^Xa) - (Xa^a) - (a^aX)) . 



(43) 



We thus find the covariance matrix of the cavity-nuclei system to be 
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(44) 



where 



_7t 



m 



n 



k 



^sinh(2z.t)+(| + — 



7 



7 



cosh (2z/t) + % 



(45 a) 



l + 32^e" 



sinh [yty 



7* 
2 



■ sinh {yty 



92 



sinh {2vt) 



(456) 
(45 c) 



According to 

Tri = STr.K-rS'^ 



' 921 

L2z/2 • ■ z/ 

there exists a symplectic transformation S (cf. Appendix A ) such that 
= diag(A^, AJ, Ag, A2) where {A^,A2} are the symplectic eigenvalues of 



This allows us to calculate the covariance matrix of the pure nuclei-cavity-output 



field system 



r 

ns— c— o 



S'Tj^iS'' 



(46) 
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gEoF 




Figure 9. Plot of the Gaussian entanglement of formation (gEoF) of the nuclei with 
the output field vs t for different values of 7/5. At 92^1 = 7 the coupling is switched 
off. The curve saturates when all excitations have leaked out of the cavity. 



with Fd' in 2 X 2 block-matrix form 

f cosh(2ri)l2 sinh(2ri)cr2 

cosh(2r2)l2 sinh (2r2)cr2 

sinh(2ri)cr^ cosh(2ri)l2 
y sinh(2r2)(T^ cosh(2r2)l2 J 



(47) 



where cosh ri = and cosh r2 = A2 and S" = | ) . One of the symplectic 

l4x4 

eigenvalues {AJ, A^} is 1, indicating a pure - and therefore unentangled - mode in the 
system. That implies that there is a single "output mode" in the out-field of the cavity to 
which the cavity-nuclear-system is entangled and we can thus trace out the unentangled 
output mode. 

The procedure for entangling the nuclei with the output field (write-in) is: let ifgq 
act for time ti to create a two-mode squeezed state 4'{g2,ti)- nuclei entangled with 
cavity and output field. To obtain a state in which the nuclei are only entangled to the 
output field, we switch the driving laser off ((72 = 0) and let the cavity decay for a time 
h 3> Tcav, obtaining an almost pure two-mode squeezed state of nuclei and the output 
mode. We define the coupling as 

9t= {g2,t<ti0,t>ti (48) 

For the parameters used in Section [3l g2 ~ 1.9 ■ 10^^/ieV. The entanglement of the 
different subsystems can be quantified: We compute the Gaussian entanglement of 
formation (gEoF) |l3] of the reduced covariance matrix of the nuclei-output field-system 
to quantify the entanglement of the nuclei with the output field (see figure [9]). The gEoF 
measures how costly it is to generate a state by mixing pure Gaussian states. It gives 
an upper bound to the entanglement of formation (EoF) and is in the present case 
equivalent to the logarithmic negativity [53j. The entanglement of the pure cavity- 
nuclei-output mode-system can be quantified using the entanglement entropy Se [Ml- 
We plot Se for the nuclei-cavity system with the output mode [see figure fTOk)] and of 
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(a) (b) 



Figure 10. a) Plot of the entanglement entropy Se of the nuclci+cavity with the 
output field vs t for different values of 7/(7. At 32^1 = 7 the coupling is switched off. 
b) Plot of the entanglement entropy Se of the nuclei with the cavity+output field vs 
(72 i for different values of 7/32- 



the nuclei with the cavity-output mode-system [see figure [TOb)]. The entanglement is 
plotted versus g2t for different ratios of the cavity decay constants and the coupling, 

7/^2- 

6.2. Write-in: Teleportation channel 

The entangled state between nuclei and the cavity output field allows us to map a state 
of a traveling light field to the nuclei using teleportation (see figure El) [26] . 

To realize the teleportation, a Bell measurement has to be performed on the output 
mode of the cavity and the signal state to be teleported. This is achieved by sending 
the two states through a 50:50 beam splitter and measuring the output quadratures 
[26] . To be able to do this, we need to know Bq, the output mode of the cavity. In the 
following, we derive an exact expression for this mode. 

We fix a time t and denote hj B{y,t), y ^ N a. complete set of bath modes outside 
the cavity. B{y,t) can be expressed as a superposition of bath operators c{x,t) 

B{y,t) = J z{y,x,t)c{x,t)dx (49) 

where we introduce a complete set of orthonormal mode functions z{y,x,t). The bath 
operators c{x, t) are known from the input-output relations [M] 



c(a;,t) = — a(t - x)x[o,t](a;), (50) 

where a{t) is given by (139 al) . To calculate B{y,t) we thus need to determine z{y,x,t). 
This can be done, calculating the variance {c^{x,t),c{x',t)) = {c^{x,t)c{x',t)) - 
(c"'^(x, t)) (c(x',t)) following two different pathways: With ( |50l) we find 

(ct(x, t), c(x', t)) = Jg(t - x')q{t - x)* , (51) 
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where q{t) is given by (1411) . Another way to express c{x,t) follows from 

c{x,t) = Y,4y,^,trB{y,t). (52) 

y 

As shown in Section \67L\ there exists only one output mode which we label y = 0. This 
mode contains all the output photons. Therefore (^B{y,ty B{y' ,t)'^ = K dyodyiQ and the 
variance using ( !52l) reads 

{c\x,t),c{x\t)) = K z{0,x,t)z{0,x',ty. (53) 

Comparing (ISTj) to (!53|) we find 

z{0,x,t) = ^M^^]L= (54) 
/ \q{t — x)\'^dx 



and K = |y J \q(t — x)\'^dx and we have thus fully determined B(0,t) (see figure [T3|) . 
Note that the bath modes are given in a frame rotating with ui + uJi to which we 
transformed in Section [2] (cu^) and Section [6] (cui). 

Therefore a state of a traveling light field can be teleported to the nuclear spins 
up to a random displacement that arises from the teleportation protocol [55l [26] . The 
random displacement can be undone, letting the beam-splitter interaction Hbs [given 
by ( fTOi) ] act for a short time, while pumping the cavity with intense coherent light as 
suggested in [50] . 

Next, we want to consider the quality of the teleportation. Whereas before (see 
figures [9] and [TOl) the time evolution of the system for a fixed switch-off time g2ti = 7 
was considered, we now consider the "final" entangled state of nuclei and output field 
depending on g2ti, where the cavity has decayed to the vacuum state while the nuclei 
are (still) stationary. 

The fidelity with which a quantum state can be teleported onto the nuclei is a 
monotonic function of the two-mode squeezing parameter 

ri = -arccosh(m(t = ti)) (55) 

with m defined in (145 al) . A typical benchmark [39] is the average fidelity with which an 
arbitrary coherent state can be mapped. This fidelity has a simple dependence on the 
two-mode squeezing parameter ri of the state used for teleportation and is given by [56] 

^ (56) 



'tel 



1 + e"2^i ■ 

We plot the teleportation fidelity dependent on the switch-off time ti (see figure [TTl) . 

Already for ri{ti) ~ 1 fidelities above 0.8 are obtained. After switching off the 
coupling we have to wait for the cavity to decay which typically happens on a nanosecond 
timescale and does not noticeably prolong the protocol. 
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Figure 11. Plot of the teleportation fidelity vs 32^1 for different values of 7/32- 



6.3. Read- out 



The beamsplitter Hamiltonian ifbs [given by ffTOj) ] enables read-out of the state of the 
nuclei by writing it to the output field of the cavity. The quantum Langevin equations 
of cavity and nuclear operators lead to almost identical solutions as for Hsq [see (139 ap 
and ( 13961) ]: of course, now a{t) is coupled to b{t) instead of b'^it) but the only other 
change to (139 ap and fl39fep is to replace u by 



'il/Ay-gl (57) 
This has the effect that all terms in (139 ap and fl395p show exponential decay with t. The 



-221; 



decay of the slowest terms ~ e~ ~ sets the timescale for read-out. To calculate the 
read-out fidelity, we need to know the state of the output field at time t = T. We 
assume that the state we want to read-out is a coherent state with displacement a^s 
at time t = fully described by its covariance matrix 76(0) = 1 and its displacement 
db{0) = {b) = «ns (while cavity and output field are in the vacuum state at t = 0). As 
the norm of the displacement ||(i(t)|| of the nuclei-cavity-output system 



(58) 



does not change under the beamsplitter transformation, the displacement of the output 
mode Bq is given by 

\dBo{t)\ = V\\d{0)r-da{t)^-d,{t )^ 

= Vi-MtW + \p4mKs\ (59) 

where q{t) and p+(t) are defined by (139 ap and fl39fep with z/ replaced by z>. 





( da{t) ] 




( m) \ 


d{t) = 


db{t) 




m) 




\ dBoit) J 




{ (Boit)) J 



At finite times, the nuclear excitations and the cavity have not fully decayed 
which leads to a loss of amplitude of the mapped state. The loss is very small for 
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Figure 12. Plot of the read-out fidelity vs git for different values of j/gi 

sufficiently large T. To assure high fidelity even for states with large photon number, 
we can amplify the output field as in Section 14.11 Then the state of the output field is 
{l Bold Bo) = («:l,ans) with n as defined in Section |4?T1 This leads to a read-out fidelity 
(see figure [12]) given by 

-^read = | ( 1, Ons | Ons ) |^ = 1 " {\q\^ + IP+H, 

where we have used relations for the transition amplitudes (as in Section I4.ip given by 

m- 

6.4- Output mode 

In figure[T3]we plot the output mode of the cavity given by fj49|) for write-in and read-out, 
respectively, and for several choices of the parameters gi^2 and 7. We are considering 
here only the idealized case of a one-sided and one-dimensional cavity. In general, the 
actual geometry of the cavity at hand has to be taken into account to determine Bq. In 
the following we briefiy discuss the shape of the mode-function. It provides some insight 
into the dynamics of the mapping process, since due to fl50|) the weight of c{x, T) in 
-8(0, t) refiects the state of the cavity mode at time t — x in the past. 

Write-in: Let us consider the two extreme cases of very strong and very weak 
cavity decay. In the former case (7 ^ 5^2) the cavity mode can be eliminated, i.e., the 
nuclear spins couple directly and with constant strength ~ g'^/'-^ to the output field: 
Zq is a stepfunction which is for g2 = ^ and constant otherwise. This is refiected in 
figure fT3l where for 7 = 100(^2 most of the excitations decay directly to the outputmode 
such that zq takes a "large" value at the time the squeezing is switched on and then 
increases only slowly in time. After switching the squeezing interaction off the cavity 
quickly decays to the vacuum. For •y <ti g2, instead, two-mode squeezing builds up in 
the nuclei-cavity system as long as the squeezing interaction is on (3yus in figure [T3|) and 
after g2 is switched off the cavity decays to its standard exponential output mode. The 
intermediate cases in figure [T3h) show the shifting weight between "initial step-function" 
and subsequent exponential decay. 
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Read-out: In the case of the beamsphtter interaction, the same cases can be 
distinguished. For large 7/(71, the cavity can be ehminated and the nuclear spins are 
mapped directly to the exponential output mode of a cavity decaying with an effective 
rate gf/'j- For smaller 7, the output mode reflects the damped free evolution of the 
nuclei-cavity system, which in this case includes oscillations (excitations are mapped 
back and forth between nuclei and cavity at rate 5^1) in absolute value and phase. 

6.5. Linear Optics with the Nuclear Spin Mode 

The interaction we have described can not only be used to map states to the nuclear 
spin ensemble but also for state generation and transformation. In fact, from a nuclear 
spin mode in the vacuum state, all single mode Gaussian states can be prepared. To 
see this, we have to show how any desired 2x2 correlation matrix F and displacement 
d E <C can be obtained. 

As we remarked already when discussing the write-in via teleportation, the beam- 
splitter Hamiltonian i^bs can be used to realize displacements of the nuclear mode. 
Driving the cavity mode with a strong laser to a coherent state with amplitude a (and 
the same phase as d) and switching on i^bs for a time t = \P\/{gi\a\) provides in the 
limit of large a a good approximation to the displacement operation by /3 |50j . 

Concerning the CM, we use that every CM of a pure Gaussian state is of the form 
F = ODO^, where D is a positive diagonal matrix with determinant one and O is 
orthogonal and symplectic. O can be seen as the effect of time evolution under some 
quadratic Hamiltonian acting on the single- mode squeezed state with CM D. In the 



Figure 13. The output mode in one dimension: Plot of zq vs position x, where x = 
is the position of the cavity a) Write-in: The squeezing interaction is "on" for 3/iS and 
then switched off. b) Read-out: For 7 ^ (?2 the excitations do not have fully decayed 
to the output mode after t = 12/is. The read-out fidelity given by (|60|) corresponds to 
the probability that the excitations in the nuclear spins have decayed into the output 
mode of the cavity. For 7/(71 = 1 and 7/(71 = 10 the read-out fidelity is -Fload > 0.98 
after t « 16 — 20/iS. For 7/51 = 100 however, it takes 200/iS to achieve F^^ad > 0.98. 
Note that for input and output modes that have similar shapes (e.g. for a network), 
it is best to consider the case where j/gi 3> 1.) 




(a) 



(b) 
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single mode case, any O represents a phase shift and is obtained by letting the nuclear 
system evolve "freely" (without laser coupling, i.e. a polarized electron interacts off- 
resonantly with the nuclei) according to the Hamiltonian oc b^b for some time. Thus the 
state with CM F can be generated in a two-step process: first generate the state with 
T = D, then apply O. 

While in the preceding paragraphs we could show how to realize operations that can 
act on any input state, no such possibility seems to exist for squeezing in our context. 
Instead we show how to obtain the pure single mode squeezed state with CM D from the 
vacuum state. Letting H^q act on the vacuum results in a two-mode squeezed state with 
squeezing parameter r2. Performing a homodyne measurement (of the X quadrature) 
on the optical part of this state projects the nuclear system into a squeezed state with 
squeezing ri = ln[cosh(2r2)]/2 [57], thus given enough two-mode squeezing, any CM D 
can be produced. 

One can go even further and simulate evolution according to any quadratic 
Hamiltonian on the nuclear-optical system: According to [HE], the Hamiltonian given by 
([8]) with the interaction part giob'' + g2a'^b'^ + h.c. enables simulation of any Hamiltonian 
quadratic in a,b,a\b^ . 

7. Remarks on internal nuclear dynamics and approximations 

With regard to the realization of the proposed protocol and the applicability of the 
approximations leading to the Hamiltonians ffTOl) . ffTTj) there are three aspects to 
consider: spontaneous emission of the quantum dot, the internal nuclear dynamics and 
errors in the bosonic description. We assume the strong coupling limit of cavity-QED 
and neglect spontaneous emission of the quantum dot. The other two aspects will be 
studied in the following. Note that the results on the internal nuclear dynamics are 
corroborated by independent work of Kurucz et al. [59]. They introduce the bosonic 
description to analyze the performance of a nuclear spin quantum memory and show 
that the performance of the memory is enhanced due to a detuning between excitations 
in the mode b versus those in other modes fe^^o and that secular dipolar terms do not 
affect the memory. 

7.1. Internal nuclear dynamics 

Up to now, we have focused exclusively on the hyperfine interaction and neglected 
"internal" nuclear dynamics, dominated by dipolar and quadrupolar interactions. 
Moreover, the hyperfine coupling leads to a dipolar interaction between nuclei mediated 
by the electron. We study the dipolar interaction between nuclear spins which is 
significantly weaker than gn, gi and g2'. the energy scale for dipolar interaction between 
two nuclei has been estimated ~ lO^^/ieV for GaAs [27J. However, since for 10"^ — 10^ 
nuclei there are many of these terms, they might play a role at the 10 — 50/is time scales 
considered. 
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7.1.1. Dipolar interaction The Hamiltonian of the direct dipolar interaction between 
nuclei is given by 

£5 E E (5(IlMl£M) _ 1,1,) , (60, 

where Yij is the vector connecting spins i and j and /Xj = is the magnetic 

moment of the nuclear spin operator Ij. Hdd can be written as 

N N 

Hm=Y.11 "i^M^^It^l + B.,ltl- + + A,/- 17 + h.c.)] (61) 

i=l j^i=l 

where Aij = 1 - 3cos^%, = -^(1 - 3cos^%), Cij = -| sin^ 6'jje~^*'^'^ Aj = 
— |sin6'ij cos6'ije*'^*-' and = /io/ij/ij/47rr?-. In GaAs the nearest-neighbor dipolar 
interaction strength is around 7 = 10~^/ieV [27]. We want to calculate the strength 
of the dipolar interaction between the main bosonic mode (that is defined as the 
mode that is coupled to the electron spin) and other bath modes (here, we no longer 
assume homogeneous coupling of the nuclei to the electron). We therefore write the 
Hamiltonian in terms of collective nuclear spin operators, use, in a next step, the bosonic 
approximation and finally separate the relevant terms (the ones which couple the main 
bosonic mode to bath modes) and calculate the coupling strength of the main mode to 
the bath modes. 

For highly polarized nuclear spins, the first term of Hdd can be written as 

N N ^ N N , K 

E E ^ 2 5^ S "'-^^'^ 2 " ^^^^ - ^^^^ ' 

i=\ j^i=l i=l j^i=l ^ ^ 

where we write If = — 1/2 + I^Ij^ and neglect the second order term I^I^I^IJ which 
requires two excitations to be non-zero; thus in the highly polarized case the contribution 
from these terms is by a factor of p = (1 — P)/2 smaller than the terms we keep. The 
last term is (for spin 1/2-nuclei) 

N N ^ N N 

E E l^MIJ ^ - 2 E E (63) 

neglecting higher order terms. In extension to the definition of the collective operators 
in Section [21 which we now label Aq, we introduce a complete set of collective 
operators Af, = J^i'^i h with = 0, .., — 1 with an orthogonal set of coefficients 
af''^ for which = 1 for every collective mode k. Defining a unitary matrix U with 

columns ck^'^^ ^ \ ^^tv^) can write 

(/r,..,/^)^ = t/A- (64) 

where A" = diag ( . ^ 2 , .., ==f I (Ag , ..,v4^_J'^. Writing Hdd in terms of 

the collective operators A'^'^ and neglecting higher order terms, 

Hdd = A+U^SUA- + {A-U^MUA- - \dUA- + h.c). (65) 
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Here, Mij = 'jijCij for i ^ j, Mij = for i = j, Sij = 'jijEij for i ^ j and 
Sii = Y.f=i^ii^ii for i = j- -D is a vector with entries Dj = Yl!ijLj=ilijDij. Next, 
we write Hdd in terms of bosonic operators, using the bosonic approximation introduced 
in Section[2l and map A~ — ^ b = (60, ••, bjy-i)'^ . This allows to separate relevant terms 
of Hdd, which couple the main bosonic mode Bq to other (bath) modes bk- Isolating the 
terms containing bo, we find 



+ h.c. 



^(f/t2M[/)ofc6, + {U^SU)okbi - \DkUok 

+ {U^SU)ooblbo + {U^2MU)oobobo + h.c. (66) 

where the notation {U'' SU)qi denotes the element (0,/) of the matrix WSU. The first 
term describes the passive coupling of the main mode 60 to other modes b^ and acquires 
a factor of two as the terms that describe the active coupling in fl66l) can be written 
bi{WMU)ikbk + bk{U^MU)Mbi = bi{W2MU)ikbk as {WMU)m = {U^ MU)iu: the entries 
of U are real so that {U^Y = [/ and M = M^, i.e., Mij = -^f sin^ Oije-^'^'^ = Mji as 
(pji = 71 + (pij. The second term in (l66ll describes the passive coupling of 60 to the modes 
bl and the third term displaces the main mode. The last two terms describe a constant 
energy shift (~ &0&0) and a squeezing term (~ 60^0 + h.c), respectively. 

The terms that couple the main mode 60 to bath modes can be written as 



60 I J](f/t2Mf/)ofc&fe + {U^SU)okbi 




rj 

2 

kj^O 

where the linear combinations of bosonic modes bk, b^ can be transformed to bosonic 



bo I C161 + C2bl - 2 > . Wok I + h.c. (67) 



modes bi and b\. The coupling strength of 60 to the first term in fl67|) is given by 
[cifei, (ci6i)t] = |ci|2 = ^ |(f/t2Mf/)o,p 

= (f/t4MMtt/)oo - |(f/t2Mf/)ooP = (AMo)^ (68) 

and |c2p = (A^o)^ for the second term. AMq and A^o depend only on the electron 
wave function and the lattice geometry. To numerically calculate AMq and ASq and 
the effect of the last two terms in fl5^ . we consider the case where the nuclei lie in a 
2-dimensional square plane with length R = VNrQ of each side on a grid with equal 
spacings tq (=0.24nm in GaAs [27])[see figure [TUfa)]. Consequently, 9ij = 7r/2, which 
simplifies many expressions in Hdd- These assumptions can be made as the height of 
the QD is small compared to its diameter, so that the variation of 6 that is dependent 
of the height of the QD is small, % ^ tt/2. 

To illustrate our results we consider two simple choices for the electron wavefunction 
such that = ^j^^^fiM^i) with = {xi,yi), 

/■(^')-"^(il)^--(il)' (-) 
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and 



exp 



(70) 



To show that the direct dipolar interaction is a weak effect compared to the optical- 
nuclear coupling g, we calculate the ratios 



di 



7(AMo) 7(AMo) 8A'j,_Ue ^ (AMq) 



92 



8A' lie 



N 



a 



(0)2 



(71) 



and d2 = • For the parameters used for the simulation in Section [3l 



L 



4 ■ 10~^ with 7 for GaAs [27]. A plot of di and d2 is shown in figure [MT b). di and d2 
are both on the order of 10~^ — 10~^, for N > 1000 nuclear spins and increase slowly 
with A^. The last two terms in (|66l) . {W SU)QQblbo and {W2MU)oobobo are small and 
zero, respectively, as can be seen in figure [T^ c): The ratio of ^3 = '^^^''^^'loo jg |-]^g 

IS zero 



order of 10 ^ — 10 ^ for A^ > 1000 nuclear spins and the ratio d^ 



92 

7([/t2Mj7)oo 

due to the symmetry of the electron wavefunction in this setting. Shifting the electron 
wavefunction such that it is not longer symmetric with respect to the coordinate origin, 
d^ is on the order of 10~^. We assume that the nuclei lie in a plane, so there is no 
displacement of 60 as Dij = for 9ij = |. Therefore, we have shown, that direct dipolar 
coupling is an effect that does not affect our protocol. 

The hyperfine coupling between electron spin and nuclear spins leads to a mediated 
dipolar interaction between nuclear spins [61j. In the bosonic description, the electron 
couples solely to the 60 mode, thus, the mediated coupling leads only to an energy shift 



^2 



(72) 
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Figure 14. a) Cosine-shaped wavefunction of the electron on a 2-dimensional 
square grid with the nuclear spins located at the vertex points, b) plot of the ratios 
'^1 — ^'''^j^^""' ^-iid d2 — ^^•'^5°'' foi' ^ cosine and a Gaussian shaped wave function. 
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For N — 10 , both ratios di and d2 are on the order of 10 , together with c), we 



see that the dipolar interaction is negligible, c) plot of the ratios ^3 



j{U^SU)oo 



and 



(i4 — ^'■'^ for a cosine and a Gaussian shaped wave function. For N — 10"*, d. 



-3 IS 



on the order of 10 ^ and di is zero due to the symmetry of the electron wavefunction. 
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that depends on the Zeeman sphtting (Dg and the number of nuclear excitations. This 
was aheady present in ([8]) and is not affecting the protocol, in fact it can help as Kurucz 
et al. [62] showed. 

For spin-1/2 systems, as considered here, the quadrupolar interaction is not present. 
For large spin I (e.g. 3/2 or 9/2) nuclei present in GaAs, there is a significant 
quadrupolar term. Depending on the strain, up to gg < 10~^/ieV have been measured 
|63j . Therefore, for / > 1/2, dots with small strain have to be considered. The 
quadrupolar interaction [60] can be treated on a similar footing as the dipolar coupling 
in Section 

7.2. Errors in the bosonic picture 

We have relied on a simple bosonic description of the collective nuclear excitations 
and neglected all corrections to that simplified picture. For homogeneous coupling 
{aj =const) this is the well-known Holstein-Primakoff approximation [33] and for 
systems cooled to a dark state [61] at moderate polarization {{A^) on the order of 
— 1/2) spin, replacing the collective spin operators by bosonic operators is accurate to 
o{l/N). The generic inhomogeneous case is discussed in detail in [35]. In that case, 
the Hamiltonian ([3]) can be seen as a zeroth order approximation in a small parameter 
~ g(l — P), where g > 1/2 and g = 1/2 for a homogeneous wave function. The first- 
order correction analyzed in [35] contains two contributions: (i) a polarization dependent 
scaling of the coupling-strength gn which has negligible effect on the adiabatic transfer 
we consider and (ii) an effective coupling of b to bath modes due to the inhomogeneity 
of the term. This correction can be computed similarly to the one in the preceding 
subsections by rewriting A^ in terms of bosonic operators. The coupling strength of the 
leading term is found to be ~ A/N = gn/VN and is thus much weaker than gi/2. Since 
5^1/2 also characterizes the energy splitting between different excitation-manifolds in the 
JC system, this term is further suppressed by energy considerations. 

8. Summary and Conclusions 

We have shown how to realize a quantum interface between the polarized nuclear spin 
ensemble in a singly charged quantum dot and a traveling optical field. The coupling 
is mediated by the electron spin and the mode of a high-Q optical cavity to which the 
quantum dot is strongly coupled. Our proposal exploits the strong hyperfine and cavity 
coupling of the electron to eliminate the electronic degree of freedom and obtain an 
effective coupling between cavity and nuclei. First, we have studied several possibilities 
to directly map the state of the cavity to the nuclei and discussed error processes and 
drawbacks of these schemes. Then, we have presented a more sophisticated interface 
which is robust to cavity decay. Read-out is achieved via cavity decay while write-in 
is based on the generation of two-mode squeezed states of nuclei and output field and 
teleportation. For typical values of hyperfine interaction and cavity lifetimes, several 
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ebit of entanglement can be generated before internal nuclear dynamics becomes non- 
negligible. All proposed schemes take advantage of the bosonic character of the nuclear 
system at high polarization, which implies that all the relevant dynamics of nuclei, 
cavity and output field is described by quadratic interactions. This allows the analytical 
solution of the dynamics and a detailed analysis of the entanglement generated. We show 
that apart from mapping a light state to the nuclei, the couplings described enable the 
preparation of arbitrary Gaussian states of the nuclear mode. 

For highly polarized nuclear spin systems the bosonic description provides a very 
convenient framework for the discussion of (dipolar and quadrupolar) "internal" nuclear 
dynamics. It is seen that these processes do not appreciably affect the performance of 
the interface. 

Our results give further evidence that nuclear spins in quantum dots can be a 
useful system for quantum information processing. In view of the recent impressive 
experimental progress in both dynamical nuclear polarization of quantum dots and 
quantum dot cavity-QED, their use for QIP protocols may not be too far off. 
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Appendix A. Gaussian states and operations 

Gaussian states and operations play a central role in quantum information with 
continuous variable systems [65]. To make this work self-contained we briefly summarize 
here the main properties of Gaussian states and operations with particular regard to 
their entanglement. 

Gaussian states are a family of states occurring very frequently in quantum optics, 
e.g., in the form of coherent, squeezed, and thermal states. Despite being defined on an 
infinite dimensional Hilbert space [JF+(]R^^), the symmetric Fock space over R^^] they 
are characterized by a finite number of real parameters, namely the first and second 
moments of pairs of canonically conjugate observables {Qi, Pi, . . . ,Qn, Pn) = R- 

One way to define them is that their characteristic function, i.e., the expectation 
values x(0 = ti'(W^cP) of the displacement operators = exp{i^^R),C^ G is a 
Gaussian function |66j : 



The displacement vector d G H and the 2N x 2N real positive covariance matrix 
(CM) 7 are given by the expectations and (co)variances of the Rk- 



X{p) = exp(-z^' 



■^d-l/Af^O- 



(A.l) 



dk = tY[pRk], 

iki = {RiRj + RjRi) — '^{Ri){Rj)- 



(1.2a) 
(1.26) 
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All d G are admissible displacement vectors and any real positive matrix 7 is a 
valid CM if it satisfies 'j > icT^ when the symplectic matrix is 

(TN = with f^l = ^ ° ^^'^^ 

The last condition summarizes all the uncertainty relations for the canonical operators 
Rj. These operators are related to the creation and annihilation operators aj, aj by the 
relations Qj = {aj + a}-)l\pl and Pj = —i{aj — aj)/v^. 

An example for a one-mode Gaussian state is a coherent state |a), with covariance 
matrix 7 = 1 and displacement d = (Re|a|,Im|Q!|)/v^. 

Entanglement: All information about the entanglement properties of Gaussian 
states is encoded in the CM. Given a CM, there are efficient criteria to decide whether 
a Gaussian state is entangled or not. 

To apply these criteria, it is useful to write the CM of a bipartite N x M Gaussian 
states in the following form. 




7 = „ , (1-4) 



where the 2A^ x 2N (2M x 2M) matrix A (B) refers to the covariances of the quadrature 
operators associated with the first (second) system and C contains the covariances 
between the two systems. A (B) are the CM of the reduced state in the first (second) 
system only. 

In the case of a two-mode system the criteria [671 168] are necessary and sufficient 
for separability: a state with CM 7 is entangled if and only if det 7 + 1 — det A — det B + 
2detC ^ 0. In this case, entanglement is necessarily accompanied by a non-positive 
partial transpose (npt) [69]. For more modes, entangled states with positive partial 
transpose exist [70] and more general criteria to decide entanglement have to be used 
[TUli]. 

For pure states, the analysis of entanglement properties becomes particularly easy 
since all such states can be transformed to a simple standard form, namely a collection of 
two- mode squeezed states (TMSS) and vacuum states, by local unitaries [73] , hence the 
entanglement of such a state is fully characterized by the vector of two-mode squeezing 
parameters. This also shows that for a 1 x M system in a pure state one can always 
identify a single mode such that only it (and not the M — 1 other modes) is entangled 
with the first system. 

For many Gaussian states it is also possible to make quantitative statements about 
the entanglement, i.e. to compute certain entanglement measures. For pure N x M 
states, the entropy of entanglement can be computed from the symplectic eigenvalues 
of the reduced CM A (or, equivalently, B). These are given by the modulus of the 
eigenvalues of <7j\fA [74j. All symplectic eigenvalue A > 1 corresponds to a TMSS 
with squeezing parameter arccosh(A)/2 in the standard form of the state at hand and 
contributes log2 A^ — (A — 1)^ log2(A^ — 1) to the entanglement entropy of the system. 
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For mixed states, it is possible to compute the negativity [71] for any N x M system 
from the symplectic eigenvalues of the CM of the partially transposed state (which is 
related to the CM obtained by replacing all momenta Pj in the second system by —Pj)- 
Every symplectic eigenvalue A < 1 contributes — log2 A to the negativity. 
For 1x1 Gaussian states with det A = det i? (so-called symmetric states), the 
entanglement of formation (EoF) can be computed [75] and for more general states 
a Gaussian version of EoF is available [13]. Even if the states are not certain to be 
Gaussian, several of the Gaussian quantities can serve as lower bounds for the actual 
amount of entanglement [76] . 

Gaussian operations: Operations that preserve the Gaussian character of the 
states they act on are called Gaussian operations [57]. Like the Gaussian states they are 
only a small family (in the set of all operations) but play a prominent role in quantum 
optics, since they comprise many of the most readily implemented state transformations 
and dynamics. With Gaussian operations and Gaussian states many of the standard 
protocols of quantum information processing such as entanglement generation, quantum 
cryptography, quantum error correction and quantum teleportation can be realized [65j . 

Of particular interest for us are the Gaussian unitaries, i.e. unitary evolutions 
generated by Hamiltonians that are at most quadratic in the creation and annihilation 
operators. Unitary displacements are generated by the linear Hamiltonian C,^R. All 
other Gaussian unitaries can be composed of three kinds [77], named according to their 
optical incarnations. The phase shifter {H = a^a) corresponds to the free evolution 
of an harmonic system. The beam splitter {H = + a^) couples two modes. Both 
generators do not change the total photon number and are therefore examples of passive 
transformations. The remaining type of Gaussian unitary is active: the (single-mode) 
squeezer is generated by the squeezing Hamiltonian H = + {a'^Y , which, when acting 
on the vacuum state decreases the variance in one quadrature {Q) by a factor / < 1 
and increases the other one by 1//. Combining these building blocks in the proper way, 
all other unitaries generated by quadratic Hamiltonians, e.g. the two-mode squeezing 
transformation [H = ab + a^b^) can be obtained. 

Both active and passive transformations map field operators to a linear combination 
of field operators (disregarding displacements caused by linear parts in the Hamiltonians, 
which can always be undone by a further displacement), i.e. for all Gaussian unitaries 
we have in the Heisenberg picture 

URU^ = SR = R'. (1.5) 

Here S' is a symplectic map on R^^, i.e. S preserves the symplectic matrix ctat, assuring 
that Ri and R[ satisfy the same commutation relations. We denote by Us the unitary 
corresponding to the symplectic transformation S. Passive operations correspond to 
symplectic transformations that are also orthogonal. 

In the Schrodinger picture. Us transforms the Gaussian state with CM 7 and 
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displacement d such that (7, d) i-^ {S'-yS'^, Sd). The two-mode squeezing transformations 

T(r)=( ^°^^^'^)^ sinh(r)a,. \ 
I sinh (r)(Ta; cosh(r)l / 

used in Sec. 14.11 is an important example of a active symplectic transformation. 

Besides Gaussian unitaries, Gaussian measurements are another important 
and readily available tool. Gaussian measurements are generalized measurements 
represented by a positive-operator- valued measure { I7, c?)(7, c?| 7, d G R^^} that is 
formed by all the projectors obtained from a pure Gaussian state I7, 0)(7, 0| 7, by 
displacements. The most important example is a limiting case of the above: the 
quadrature measurements (von Neumann measurements which project on the (improper, 
infinitely squeezed) eigenstates of, e.g., Q). In quantum optics, these are well 
approximated by homodyne detection. For example, the "Bell- or "EPR-measurement" 
that is part of the teleportation protocol is a measurement of the commuting quadrature 
operators Qi + Q2 and Pi — P2- 



Appendix B. Landau-Zener transitions 

In a rotating frame with U = exp[—^ f{uJi — uj2)cTz + {uji + U2)ldt] the Heisenberg 
equations, given by ( fT2l) . read: 



u' = — igexp yi J {uji — uj2)dtjv' (2.1) 
v' = — igexp I —i / {ui — U2)dt \ u' (2.2) 



The initial boundary conditions of the coupled differential Equations (12.11) and (12.21) are 
now chosen such, that the photon operator a at time t —00 is mapped to the nuclear 
spin operator b at t —>■ (yo 

u'_^ = l,\v'_J=0. (2.3) 

Eliminating u' in (12.11) and (12. 2p leads to the single equation: 

v' + iptv' + gh' = 0, (2.4) 

where g = 0. Together with the substitution v' = ■l^'^^^'^^^'^^Ui, (12.41) reduces to the 
so called Weber equation: 

Ul+[9'-^^ + ^t'^U, = 0. (2.5) 

Solving (12. 5p as proposed by Landau and Zener and considering the asymptotic behavior 
of the solution at t — »• 00, it is found to be 



hm U,{t) = -^ p, I,, e-^^-e-^*^(v^t)-^-, (2.6) 

t^oo 1 [27^ + Ij 

2 

where 72 = ^ and the constant K = y/^exp (— • The probability that the photonic 
operator a is mapped to the collective nuclear spin operator b is given by (ITBi) . 
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